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Abstract. The 3-body crystallization diagrams of 
C/O /Ne ionic mixtures characteristic of white dwarf inte- 
riors are examined within the framework of the density- 
functional theory of freezing. The crystallization process 
is described more accurately than in former calculations 
where the three-component system was treated as an effec- 
tive two-component mixture (Segretain et al. 1994). The 
distillation process due to neon-crystallization is found to 
occur only for the late stages of crystallization. At the 
beginning, the presence of neon plays only a minor role 
and the phase diagram resembles a pure carbon-oxygen 
diagram. The hnal phase diagram is found to exhibit an 
azeotropic point with a neon concentration xn c = 0.22, 
a carbon concentration ip = 0.78 and an oxygen concen- 
tration xo = 0, so that during the distillation process, the 
fluid crystallizes into a pure neon-carbon solid. The crit- 
ical temperature is Ta = 0.85 Tc, where Tc is the pure 
carbon crystallization temperature. We use this accurate 
phase diagram to calculate the total gravitational energy 
released during white dwarf crystallization and the related 
time delay. The final result yields At ss 2.6 Gyr, among 
which about 20% are due to the neon-distillation process. 
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1. Introduction 

Following the pionneer work of Stevenson (1980), Garcia- 



Berro et al. (1988) have demonstrated the importance of 



the shape of the crystallization diagram on the cooling his- 
tory of white dwarfs (WD). In the eighties, all the work 
has focussed on the crystallization of a pure carbon-oxygen 
mixt ure (Stevenson 1980|; M ochkovitch 1983 ; Barrat et al. 



1988 ; Ichimaru et al. 1988 ). More recently, Isern et al. 
( 1991 ) have pointed out the importance of minor elements, 
in particular neon, on the general cooling process. In order 



to examine in details the effect of these d iffere nt crystal- 
lization processes, Segretain & Chabrier ( 1993 ) have cal- 
culated the crystallization diagrams of dense binary ionic 
mixtures with arbitrary charge ratios under thermody- 
namic conditions characteristic of dense stellar plasmas. 
The diagram was found to evolve from a spindle form to an 
azeotropic, and eventually an eutectic form, with increas- 
ing charge ratio Zi jZ\ . These binary diagrams have been 
used by Segretain et al. ( 1994 ) to compute the cooling 
of a C/O/Ne/Fe WD. The crystallization processes were 
shown to bear major consequences on the cooling time and 
the luminosity function of old WDs (Segretain et al. 1994; 
Hernanz et al. 1994). In these calculations, and given the 
complexity of the calculation of multi-component phase 
diagrams, the ternary C/O/Ne or C/O/Fe mixtures were 
assumed to behave as effective binary mixtures composed 
of neon (resp. iron) and an element of average charge 
(Z) — 7. The calculations showed that neon plays a major 
role in the cooling of WD, because of a distillation pro- 
cess bound to the presence of an azeotropic point in the 
phase diagram. The resulting sedimentation of neon was 
found to have an effect on the cooling time comparable 
to the one due to the crystallization of the carbon-oxygen 
mixture itself. In this paper we have e xtended the calcu- 
lations of Segretain & Chabrier ( 1993 ) to the calculation 
of a true ternary C/O/Ne diagram. The diagram is calcu- 
lated within the framework of modern density functional 
theory (DFT) of freezing, and the plasma correlation func- 
tions are calculated within the mean spherical approxima- 
tion (MSA) (see Segretain & Chabrier 1993 for details). In 
Sect. 3, we describe the resulting diagram and we compare 
it with the previous effective two-body diagram. In Sect. 
4, we use this ternary phase diagram to derive the global 
effect of neon crystallization on the cooling of a WD. 



2. Method 

The crystallization diagram has been computed within the 
framework of the density functional theory (DFT) of freez- 
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ing. One of the main advantages of the DFT is that we 
can calculate directly the free energy difference between 
the solid and the fluid phases. We consider a mixture of 
Ni ions of charge Zi (i — 1,2,3) in a volume V at a 
temperature T. The number densities and molar fractions 
read respectively pi — Ni/V and xi = N%/N (i=l,2,3) 
with N = ^ZjiVj. The electronic and ionic coupling pa- 
rameters are defined respectively as T e — e 2 /a e kBT and 
T, = zf /3 T e , where a e = (4tt /3N e /V)~ 1/3 is the mean 
inter-electronic distance and < Z >= J^i 

Xi Zi is the av- 
erage charge. The free energy of the liquid phase reads: 



3. C/O/Ne phase diagram 

Since the pressure of the degenerate electron gas is consid- 
erably larger than the ionic pressure under the conditions 
of interest (degenerate stars interiors), fluid-solid phase 
equilibrium at constant pressure and temperature is equiv- 
alent to phase coexistence at constant electron density p e 
and temperature T, i.e. at constant T e . For fixed values of 
T e and T, the concentrations of coexistent fluid and disor- 
dered solid phases are determined by usual double-tangent 
constructions. 
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where the first term on the r.h.s. is the ideal ionic con- 
tribution whereas the second term denotes the correlation 
terms. This term is well approximated by a linear interpo- 
lation of the free energies of the pure phases (designated by 
the superscript OCP -One Component Plasma- ), for which 
analytical fits are available (Stringfellow et al. |l990|) . The 
free energy of the solid writes: 

F S = F L + AF 

AF = AF ld + AF ex 

with 
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Fig. 1. C/O/Ne crystallization diagram for T e > 7. The coor- 
dinates of the azeotropic point are T a = 0.85 Tq and xq = 0.78, 
xo = 0.0, XNe = 0.22, where Tc is the crystallization temper- 
ature of pure carbon. The bold solid lines visualize the liquidus 
and the bold dotted and dashed lines the solidus. 
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where ^2' G represents the sum over all the reciproqual 
lattice vectors (RLV) for G ^ 0. The p v denote the num- 
ber density of the v component either in the solid (S) 
or in the liquid phase (L). For each concentration x v at a 
given temperature T, AF is minimized with respect to the 
width of the gaussians a v in order to calculate Fg. The 
C M „(G) are the Fourier transforms of the direct correlation 
functions of the liquid near crystallization. These latter 
are computed within the framework of the Mean Spher- 
i cal A pproximation (MSA). As in Segretain & Chabrier 
(|1993|) , we use a reduction of the second RLV in C M „(G) 
in order to obtain a bec crystal for the one component 
syste m at T e Z 5 / 3 = 178 (Slattery et al. |l982[ Stringfellow 
et al.|l990[). 



Figure 1 shows the resulting 3-component phase dia- 
gram for the C/O/Ne mixture under thermodynamic con- 
ditions characteristic of WD interiors, i.e. 10 10 < P < 10 16 
Mbar, 10 3 < p < 10 6 gem" 3 and for T < 5.10 7 K, which 
yields T e > 7. The 3-D picture of the C/O/Ne diagram 
can be first visualized with the 2-D projections on the 
C/O, C/Ne and O/Ne planes. Note that, for sake of clarity, 
only the lower part of the general diagram (T/Tc < 1.2, 
where Tg is the crystallization temperature of pure car- 
bon), is shown. The pure carbon and oxygen mixture 
(Xmb — 0) limit crystallize into a spindle diagram, as 
expected from the small charge ratio. So does the bare 
O/Ne mixture (Xc = limit), with a larger crystalliza- 
tion temperature, since T/y e > Tg and To > Tg where 
Tj , i = C, O, Ne, denotes the crystallization temperature 
of pure carbon, oxygen and neon respectively. The pure 
C/Ne (Xq = 0) diagram exhibits an azeotropic behav- 
ior with an azeotropic point for xc — 0.78, xn e = 0.22 
and T a — 0.85 Tc- The true 3-body diagram is visual- 
ized by the bold lines, where the solid line display the 
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liquidus and the dotted and dashed lines represent the 
solidus. As shown, the azcotropic point of the three-body 
diagram is the azeotropic point of the two-body C/Ne mix- 
ture (Xo = 0), which is the lowest stable point for the 
liquid in the three-component C/O/Ne mixture. 

We also note that, for small concentrations of neon, 
and for temperatures well above the azeotropic tem- 
perature, the liquid and the solid have the same Ne- 
concentration. The mixture then crystallizes like an equiv- 
alent Ne-enriched C/O mixture. In other words, at high 
temperature, and for low concentrations of neon, the crys- 
tallization diagram is not affected by the presence of neon. 
As the temperature approaches the azeotropic tempera- 
ture, neon starts to modify the crystallization process. 
The d istillation process, as described in Segretain et al. 
( 1994 ), begins. The liquid becomes enriched in neon and 
the diagram evolves into the azeotropic shape. 

This behavior is different from the one obtained with 
the effective 2-body mixture used in Segretain et al. 
( |1994D , where the Ne-enriched process was found to start 
as soon as crystallization set in. As described above, the 
present complete 3-body calculations show that C/O crys- 
tallization first takes places without Ne-enrichment. 

This bears important consequences for the astrophysi- 
cal application, as will be shown in the next section, where 
the diagram will be examined for C/O/Ne concentrations 
characteristic of WD interiors. 



4. Consequences on white dwarf cooling 

The mass fractions characteristic of WD interior are 
X Ne = 1.0% and X c = X = 49.5%, uniformly dis- 
tributed in the liquid phase (the effect of an initial non- 
uniform profile has been examined in details in Segretain 
at al. 1994). The exact C/O mass fraction in WD inte- 
riors is subject to important uncertainties. The present 
cquimolar distribution has been chosen for the sake of sim- 
plicity. Larger initial oxygen concentrations in the core of 
the WD will, at any rate, lead to even smaller energy re- 
lease by Neon sedimentation. As shown on the diagram of 
Fig. 1, the process of crystallization develops as follows. 
Crystallization begins at the center of the star where the 
density, and thus the crystallization temperature, is the 
largest (T c oc p 1 / 3 /!!%). In that mentioned in the 

previous section, we get the same behaviour as for a pure 
C/O mixture, but with concentrations Xc = Xo = 0.495, 
i.e the same neon concentration in the solid and the liquid 
phase. Until T/T c = 1.12 (point A in Fig. 2) the result- 
ing sedimentation of oxygen is the same as in a pure C/O 
WD. Beyond T/Tq = 1.12, the distillation process begins, 
so that neon concentration becomes different in the two 
phases. However, the neon effect on the release of gravi- 
tational energy will be small because, at point A, about 
70% (in mass) of the star has already crystallized. This is 
very different from the result obtained with the effective 
2-body diagram, where neon was found to crystallize first. 



At the end of the distillation process, i.e. at the azeotropic 
temperature (point B in Fig. 2), all neon has crystallized; 
the remaining liquid contains only carbon and oxygen and 
sedimentation of oxygen continues normally. The oxygen 
and neon profiles in solid along the crystallization process 
are shown on Fig 2. 
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Fig. 2. Concentration profiles at the end of the crystalliza- 
tion, for two different WDs. Dot-dashed line : oxygen pro- 
file for a pure C/O WD. Solid line : oxygen profile for a 
C/O/Ne WD. Dashed line : neon profile for a C/O/Ne WD. 
The crystallization front evolves as indicated by the arrow. 
The neon-distillation process begins at point A and ends at 
point B. During this process, the solid which forms does not 
contain oxygen, since the azeotropic point is characterized by 
Xo = 0. After point B, all neon has crystallized; the remaining 
solid contains only carbon and oxygen, and the usual oxygen 
sedimentation process proceeds again. 



These concentration profiles are used to compute the 
evolution of the binding energy along crystallization, as 
in Segretain et al. (1994). The difference of binding en- 
ergy due to crystallization is AB = 3.10 10 46 erg, only 
~20% more than the one obtained for the crystallization 
of a pure C/O WD. The effect of neon crystallization, 
when properly taken into account, is then smaller than 
when treated within the effective 2-body approximation. 
The release of binding energy due to the crystallization of 
the C/O /Ne mixture yields a total time delay about 20% 
larger than the one due to the crystallization of a pure 
C/O mixture, i.e At rs 2.6 Gyr. The effect intrinsically 
due to neon crystallization is about 0.4 Gyr. This effect 
on the cooling time, however, does not modify the age of 
the Galactic disk, based on the calculation of the WD lu- 



minosity function, obtained by Hernanz et al. ( 1994 ), as 
explain in section 3.2 of the paper of Hernanz et al. ( 1994 ). 
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5. Conclusion 

In this paper we have computed for the first time the crys- 
tallization diagram of a three- component C/O/Ne mixture 
characteristic of WD interiors. The general phase diagram 
is found to be of the azeotropic form. We use this phase 
diagram to compute the crystallization process of a uni- 
form distibution C/O/Nc WD with X Ne = 1.0% and 
Xc = Xo = 49.5%. The presence of neon is found to 
be inconsequential at the beginning of the crystallization 
process. The neon concentration remains the same in the 
solid and in the liquid phase, so that the major process 
is the sedimentation of oxygen, characterized by a spin- 
dle diagram, as found for a pure C/O WD. When 70% 
(in mass) of the WD has crystallized, the shape of the 
diagram evolves substantially, and the influence of neon 
becomes important. A distillation process begins and all 
the neon left in the liquid crystallizes at the azeotropic 
concentration/temperature, i.e T a — 0.85 X"c 7 ic = 0.78, 
xo — 0.0 and x^ e = 0.22. This crystallization process 
is different from the one obtained with the effective 2- 
body diagram for which the distillation process is found 
to begin as soon as crystallization sets in. The time delay 
produced by the afore-mentioned general crystallization 
process is about 2.6 Gyr, i.e about 20% larger than the 
one obtained for a pure C/O mixture. The effect of neon 
is only of 0.4 Gyr. These detailed calculations show that, 
though not completely negligible, the effect of neon, and 
iron as well, does not bear important consequences on the 
general cooling history of WDs. 
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